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Abstract 

We present results from our analysis of the finite-temperature properties of the spin 1/2 J1-J2 Heisenberg model on 
a square lattice. The analysis is based on the exact diagonalization of small clusters with 16 and 20 sites utilizing 
the finite-temperature Lanczos method. In particular, we focus on the temperature dependence of the third-order 
magnetic susceptibility as a method to resolve the ambiguity of exchange constants. We discuss the entire range 
of the frustration angle = tan~^(J2/Ji) parameterizing the different possible phases of the model, including the 
large region in the phase diagram with at least one ferromagnetic exchange constant. 
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The spin-1/2 Heisenberg model on a two-dimensional 
square lattice with next-nearest neighbour interac- 
tion belongs to the most intensely studied models for 
frustrated spin systems. Its Hamiltonian is of the form 



7? = Ji ^ 5, • 5, + J2 ^ 5, • Sk, 



(1) 



(«j>i 



(»'=>2 



where the sum on (ij)i runs over nearest neighbour 
and the sum (1^)2 over diagonal next-nearest neigh- 
bour bonds. We allow the exchange constants ,J\ and 
J2 to be negative (FM) as well as positive (AF) . Only 
the relative size of the exchange couplings determines 
the physics of the model. Therefore it is convenient to 
introduce an overall energy scale Jc = \/ Jl + J2 and 
a frustration angle <j) = tan~^(J2/Ji) to characterize 
the model. The schematic phase diagram of the model, 
which is shown in the left part of Fig. 1, can roughly 
be subdivided into three ordered phases, plus two spin- 
gapped phases or families of phases. The former are 
characterised by anomalies in the susceptibility x(q) 
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at the ordering vectors q = (0, 0) (FM), q = (tt, 0), or 
(0, tt) (coUinear antiferromagnet, CAF), and q = (tt, tt) 
(Neel-type antiferromagnet, NAF). 

Experimental realisations of the frustrated square- 
lattice antiferromagnet include the quasi-two-dimen- 
sional compounds Li2VO(Si,Ge)04 and Pb2VO(P04)2 Q- 
The latter is believed to have one ferromagnetic ex- 
change constant. In the following, we shall refer to this 
compound when discussing the experimental situation. 

In principle, the average interaction constant Jc can 
be determined from the asymptotic behaviour of the 
heat capacity and the magnetic susceptibility at high 
temperatures. However, the determination of the cor- 
rect frustration angle, and therefore the correct val- 
ues of the exchange constants of the compounds men- 
tioned above has proved difficult. Diffuse neutron scat- 
tering would provide a means to determine J\ and J2 
separately, therefore we have calculated the static spin 
structure factor ^(q, T) previously |2| . Combining the 
parameter dependence of the third-order susceptibil- 
ity presented here with results obtained earlier for the 
heat capacity and the linear susceptibility, our method 
provides a way to unambiguously determine the exper- 
imental values of Ji and J2 from the thermodynamic 
properties of the compounds. 
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Fig. 1. Left: Schematic phase diagram of the frustrated Heisen- 
berg antiferromagnet. The shaded areas denote the spin-liquid 
region for Ji > 0, J2 ~ Ji/'^ and a possible new spin-liquid 
region for Ji < 0, J2 ~ — Ji/2. 

Right: Temperature dependence of the third-order magnetic 
susceptibility of the frustrated Heisenberg antiferromagnet on 
a 20-site cluster obtained from exact diagonalization. The data 
are taken at two different frustration angles (p — —0.11 tt (solid 
line) and (p — 0.64 tt (dashed line). The inset shows the tiny 
minimum occurring at T S:; 2Jc/A:b- 

The third-order susceptibihty x"'{T) is defined via 
the small- field expansion of the magnetisation M{T)\ 
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where the symbol {. . .) denotes the trace over the sta- 
tistical operator for zero magnetic field B. S^^t is the z 
component of the total spin of the system, A'^ the sys- 
tem size. As usual, A'a is the Avogadro number, /xo the 
magnetic permeability, g the gyromagnetic ratio, /xb 
the Bohr magneton, and fce the Boltzmann constant. 
An expansion of the Brillouin function yields the 
high-temperature asymptotic behaviour of x"'{T), 
which is of the form 
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Note that x" {T) is negative for temperatures T ^ Jc. 
We have applied the finite-temperature Lanczos 
method [3| for clusters with 16 and 20 sites to evaluate 
the thermodynamic traces in Eq. 3. To give an exam- 
Table 1 

Characteristic values for the third-order susceptibility ob- 
tained for the two different frustration angles — — O.IItt 
and (j) — 0.64 TT which are equally possible for the compound 
Pb2VO(P04)2- The second column lists the ratio of the Curie- 
Weiss temperature to the maximum temperature of the linear 
susceptibility as obtained from Ref. ^. 
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-0.11 
0.64 


0.49 
0.49 


0.11 
0.21 


1.69 
1.45 


2.37 
2.02 


1.41 

0.64 


-0.0018 
-0.0029 



Fig. 2. Maximum value (left) and maximum temperature 
(right) of the third-order magnetic susceptibility of the frus- 
trated Heisenberg antiferromagnet. The filled dots represent 
the exact-diagonalization results for a 16-site cluster, the open 
circles correspond to a 20-site cluster. 

pie, the right part of Fig. 1 shows the temperature 
dependence of x" of a- 20-site cluster for two differ- 
ent values (f)± of the frustration angle. These values 
are chosen such that they correspond to those deter- 
mined for Pb2VO(P04)2 [EQ- The frustration angle 
(j) = —0.11 TT (solid line) corresponds to the Neel anti- 
ferromagnet, while (j) = 0.64 TT (dashed line) describes 
the coUinear phase. 

In all non-ferromagnetic phases of the model, the 
temperature dependence of x"'{T) has a pronounced 
maximum at a temperature Tmax, vanishes at a tem- 
perature Tq", passes through a tiny minimum at a tem- 
perature T^'in, and eventually approaches the high- 
temperature T~ dependence. We have followed these 
characteristic temperatures as a function of the frus- 
tration angle. The initial maximum temperature, to- 
gether with the value of x'" at that point, are shown 
for the 16-site (solid dots) and 20-site cluster (open 
circles) in Fig. 2. The maxima occur at temperatures 
Tmax ^ Jc, therefore finite-size effects are large. 

The characteristic temperatures, together with the 
values of x'" at maximum and minimum, are compiled 
in Table 1 for — —0.11 n and 0.64 tt. Also shown is the 
ratio Ocw/Tmax of the Curie-Weiss temperature to the 
position of the characteristic maximum of the linear 
susceptibility x{T) taken from Ref. |2|, which is equal 
for both values — (f)± . Together with the results from 
Ref. [3, our findings should be useful in determining 
the precise value of (f> for a given J1-J2 compound, as 
exemplified here for Pb2VO(P04)2. 
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